Sensitive predators and endurant preys in an ecosystem driven by correlated noises 
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We investigate the Volterra ecosystem driven by correlated noises. The competition of the preda- 
tors induces an increasing in population density of the predators. The competition of the preys, 
however, leads the predators to decay. The predators may have better stability under strong cor- 
related noises. The predators undergo a sensitivity to a random environment, whereas the preys 
exhibit a surprising endurance to this stochasticity. 
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Predator-prey ecosystems are one of the most amazing 
models in interacting population biology [1-3]. The prey- 
predator model and its derived ones are applied in a wide 
range of fields, for example, tumors, measles, epidemics 
and food chains [4-6] . Though predator-prey model sug- 
gests that the simple interactions can result in periodic 
behavior of the population, it is still an unrealistic as- 
sumption. To fill this deficiency, scientists struggle to 
find more realistic approaches. Stochastic models are 
so far better selections, including individual level model 
(ILM)[7] and two-species stochastic population model 
(SPM) [8-10]. The anterior probes a situation when the 
number of individuals is finite, and the later describes a 
strictly infinite one. 

Surprising phenomena emerge sometimes just because 
there exists correlation or collaboration; one of the most 
well known example is the self-organized behavior [11]. 
The population in ecosystems, including prey and preda- 
tors, are affected by simple stochastic processes of moral- 
ity, reproduction, and predation. Differential equations 
with stochastic components match their realities better 
than the deterministic ones. Due to the difficulty and 
complexity in mathematical analysis, previous researches 
in prey-predator models only consider independent noises 
[8-10]. Regrettably, the correlated noises, which are more 
realistic tentation to depict true nature than indepen- 
dent ones, are easily left in the basket. We all know that 
stochastic processes in ecosystems may come from the 
same origin like disaster, thus it is reasonable to consider 
the correlation between noises. In this paper, we focus on 
the two-species model with correlated stochastic compo- 
nents for the first time. The correlated noises we consider 
here, which distinguish from independent ones Cai and 
Lin have concerned with in their latest paper [8], induce 
some novel phenomena in a prey-predator ecosystem that 
are not found before. 

We study two-species models to elucidate the mecha- 
nism of interesting influence of correlated noises on the 
stochastic dynamical ecosystems. Consider Volterra's 
model with prey resource limitation, described by the 
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where X > and Y > denote the densities of indi- 
viduals in a population density of preys and predators, 
respectively. In the absence of noise, they will evolve 
into a stable value {Xq,Yo) as i — + oo. ri,r2 are the cor- 
responding reproduction rate. K is the capacity of prey 
population, that is, the population size for t — > cxd if 
Y = 0. The coefficient ki quantifies the impact which an 
individual predator has on the reproduction rate of an in- 
dividual prey. Conversely, k2 gives the impact which an 
individual prey has on the reproduction rate of a preda- 
tor. ^{t),ri{t) are correlated, Gaussian white noises sat- 
isfying 

mm) = 2AhS{t-t'), (3) 

mvit')) = 2M2 5it-t' ), (4) 
imvit')) - 2X^MiM2S{t-t'), (5) 

in which Mi,M2 are the intensities of noises; A, ranges 
from zero to one, denotes the correlation coefficient be- 
tween and r]{t), 6{t~t') is Dirac delta function under 
different moments. 

An approach to describing stochastic noise is using 
three important ideas: samples, events, and probability 
[12]. For example, consider a large population, say N, of 
preys. These are the samples. The events we consider are 
the collection of all possible aggregations of these preys. 
The probability of an event is simple (1/iV) times the 
number of samples in the event. Therefore, our nota- 
tion for the stationary probability distribution (SPD) of 
an event X (a population density of the preys) or Y (a 
population density of the predators) is 
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, , Number of predators in Y 
P{Y) = , (7) 

in which Ni,N2 are all preys and predators, re- 
spectively. Consider P{X) and P{Y) are inde- 
pendent, the joint stationary probability distribution 
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FIG. 1: Joint stationary probability distributions under differ- 
ent noise intensities. The parameters are Mi = 0.05, A = 0.0, 
(o) - (d) M2 = 0.1, 0.5, 2.0, 5.0. 



(JSPD) of (X,y) is P{X,Y) = P{X)P(Y) satisfying 

/o°°/„~p(xy)dXrfy = i. 

In general, the probability distribution is still a quali- 
tative parameter, though noise-induced transition can be 
depicted by observing its changing trend. Here we define 
a quantitative parameter, i.e. the mean density of pop- 
ulation, which is useful to display the influence of noises 
on the prey-predator ecosystems. The mean densities of 
the preys and predators are respectively written as 

{X) = / XP{X,Y)dXdY, (8) 
Jo Jo 

f'OC f'OQ 

{Y) = / YP{X,Y)dXdY, (9) 
Jo Jo 

and their variances are a% = (X^) — (X)^, and ay = 
(F^) - (Y)^, respectively. 

One of the approaches to solve Eq.(l) and (2) is de- 
riving the equivalent stationary Fokker-Planck equation 
[8], but that is possible only toward two non-correlated 
noises. Once consider the correlation of noises, it is diffi- 
cult to derive and solve the Fokker-Planck equation. Here 
we apply a real-time simulation of Eq.(l) and (2) to ob- 
tain the joint stationary probability distributions of the 
preys and predators densities. JSPD is useful to repre- 
sent the influence of noises on the preys and predators. 
The quantitative parameters are also calculated for deep 
analysis. For example, the mean density of the predators 
is used to display their long-term behavior. 

Figure 1 shows the peak of the joint stationary proba- 
bility distributions drops and shifts to position (0, 0) with 
M2 increasing. With an increase in Mi, similar effects 
are observed, and they are not plotted here for simplic- 
ity. Obviously, lower peak height and more large preys 
and predators densities mean less system stability. 



FIG. 2: Relationship of mean densities of the preys and preda- 
tors with the noises intensities as A = 0.0. (a) Mi = 0.05, (6) 
M2 = 0.1. 



It is comprehensible that noises can induce un-stabihty. 
Consider a long term, whether the noises will do good to 
the preys and predators or not is still not known. So 
quantitative parameters, e.g. the mean densities of the 
preys and predators are helpful to do that. Only accord- 
ing to the JSPD, it looks as if there are no differences 
between the influence of noises on the preys and those 
on the predators. Figure 2, however, illustrates their dis- 
tinctions. The average predator density increases with 
M2, and decreases with Mi. Obviously, r]{t) does good 
to the predators, but ^{t) does not. Provided that the 
former is regarded as the predators competition and the 
latter as the preys competition, it is not difficult to com- 
prehend why predators benefit from their competitions, 
and why not from the prey competitions. Surprisingly, 
the preys densities properly maintain their intactness in 
spite of the changing in Mi or M2. Why do they hap- 
pen like that? Here we suggested an important reason, 
which is easily understood from Eq.(l) and (2): both 
of the growth rate and source limitation keep the prey's 
stability, whereas the only death rate can not keeps the 
predator's stability. Perhaps this is why the predators are 
more sensitive to the stochastic noises than the preys. 

The counterpart distinctions induced by the correla- 
tion are shown in Fig. 3 and Fig. 4. Figure 3 displays 
that the stationary probability distributions of the prey 
does not change with the noise correlation. The station- 
ary probability distributions have the same fitting curve 
under different correlations. Their variances are also in- 
variable (see the inset of Fig. 3). This indicates that the 
correlation has little effect on the preys. Unlike what 
have happened on the preys, figure 4 shows that the 
peak height of the stationary probability distributions 
of the predators increases with the correlation (see the 
corresponding fitting curves), and its variance, shown in 
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the inset of Fig. 4, decreases with the correlation. Like 
what we have mentioned above, ^(i) and ri{t) refer to 
the preys and predators competitions, respectively. We 
suggest that the correlation between noises refers to the 
synchronous components of the competitions, and thus it 
is possible for the predators to maintain more stable SPD 
provided that they adjust their competition to those of 
the preys. 



FIG. 3: Stationary probability distributions of the preys un- 
der different correlations. The parameters are Mi = 0.05, 
M2 = 0.1. Inset: Variances of SPD against the correlations, 
A. 



In summary, the correlation or collaboration, which ex- 
ists in numerous dynamical systems, can cause interest- 
ing responses of the prey-predator ecosystems to noises. 
Due to the support of their growth and resource limita- 
tion, the preys undergo perfect endurance to the external 
noises. Conversely, the predators are not only sensitive 
to the noise intensity, but also impressible to the noise 
correlation. Strong correlated noises can lead the station- 
ary probability distributions of the predators to stability. 
Our model is expected to support a fact: the complexi- 
ties of biosystems originate from their collaborations and 
correlations. 



FIG. 4: Stationary probability distributions of the predator 
under different correlations. The parameters are Mi = 0.05, 
M2 = 0.1. Inset: Variances of SPD against the correlations, 
A. 
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